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Many present day quantum optics experiments, particularly in optical quantum information pro- 
cessing, rely on number-resolving photo-detection as a basic building block. In this paper we demon- 
strate that a simple classical optimisation technique can sometimes be employed to post-process the 
detector signature and improve the confidence of the measurement outcome in the presence of 
photon-number errors such as loss or dark-counts. While the regime in which this technique is appli- 
cable is rather restrictive, and will likely not be very useful for the large-scale quantum information 
processing applications of the future, the ideas presented might be employed in some present-day 
experiments where photo-detectors are typically very poor. 
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INTRODUCTION 

Number-resolving photo-detection forms the basis of 
many quantum optics experiments, particularly in the 
emerging field of optical quantum information processing 
[IHl] . Ideally we aim to employ detectors that can distin- 
guish, with high fidelity, between different Fock states. 
However, in practise, unavoidable effects such as loss and 
dark-counts corrupt the photo-detection process, lower- 
ing the fidelity of our measurement results. 

In this paper we present a straightforward technique, 
which in some parameter regimes, allows classical post- 
processing of the measurement signal to improve the fi- 
delity of the measurement outcome. We present two ex- 
ample applications for this technique, including a simple 
application of practical interest - heralded Fock state 
preparation using parametric down-conversion (PDC). 
We show that in some circumstances our technique can 
indeed improve the fidelity of the heralded state. Re- 
grettably, improvement is only possible in the regime 
where the detector is already behaving badly (e.g. lossy or 
with significant dark-counts). However, given that many 
present-day experiments employ poor photo-detection, 
our technique may nonetheless be applicable. 

We contrast our technique to the matrix inversion tech- 
niques that are used to reconstruct photon-number distri- 
butions 1^ . In such techniques we take a large number 
of measurements to determine a vector iTmeas compris- 
ing probabilities of detecting different photon numbers. 
The relationship between the actual photon-number dis- 
tribution and the measured one can be related by a ma- 
trix M, which captures information about the number- 
resolving power of the detector, including losses and dark- 
counts, iTmeas = ^'^act- Then, applying an inversion and 
fitting technique, the actual photon-number distribution 
can be inferred from a measured distribution. On the 
other hand, the technique we present here is intended 
for single shot measurements, i.e. how do we optimise a 
single measurement rather than an ensemble of measure- 
ments? This important difference makes our technique 



more applicable to quantum information processing ap- 
plications where a single measurement is used either for 
feed-forward or the measurement of the output state. 

A similar technique was discussed in by O' Sullivan et 
al. [7 , where they focussed on the ability of lossy de- 
tectors to generate Fock states via heralded PDC. While 
we discuss this example also, we differentiate the present 
work by specifically focussing on the detector optimisa- 
tion process rather than the state preparation procedure. 
In particular, we derive optimisation maps^ which show, 
for a given detector and known number distribution, how 
to map raw detector signatures to optimised detector sig- 
natures. In other words, the optimisation maps tell us 
the extent to which each measurement signature must be 
compensated. This technique may subsequently be used 
in applications broader than just heralded state prepara- 
tion. We also focus on the measurement fidelities, which 
tell us the confidence that our optimised signatures are 
correct. Finally, we demonstrate that the form of the 
known number distribution has a strong effect on the op- 
timisation map and measurement fidelities, effectively bi- 
asing the regimes in which optimisation is effective. Thus, 
the applicability of such a technique is highly application- 
dependent. 



DETECTOR MODEL 

In terms of the photon-number degree of freedom (i.e. 
ignoring spatio-spectral effects [8 ), the detector can be 
fully characterised using a matrix of conditional proba- 
bilities P(m|n), the conditional probability that m pho- 
tons were detected given that n were incident [5] fO] . This 
gives us full information about how the detector reacts 
to different input states. 

We will focus on a model consisting of an ideal number- 
resolving photo-detector preceded by a loss channel and 
subject to dark-counts. The loss rate is pioss and the dark- 
count rate ^00(^5 where i is the number of dark-counts 



[3]. Then, we can write, 
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and we assume that pdc is a Poissonian distribution [5] , 






(2) 



where A is the dark-count rate and is equal to the average 
number of dark-counts |5j. Thus, we are summing over 
ah combinations where the sum of the number of dark- 
counts, d, and the number of photons minus those lost, 
m — d, is equal to the number of measured photons m. 

In Fig. [l] we show the P{m\n) matrix for our detec- 
tor in the presence of different loss and dark-count rates. 
For a perfect detector (i.e. pioss = Pbc = 0), the matrix 
is the identity matrix since every photon number input 
state perfectly maps to the desired signature. However, 
as losses and dark-counts are introduced, off-diagonal el- 
ements emerge, corresponding to incorrect measurement 
signatures. In the case of very high loss rates the peaks 
are localised around the m = axis, i.e. we always mea- 
sure zero photons since everything is lost. On the other 
hand, in the case of high dark-count rates the peaks be- 
come very delocalised since a large number of different 
measurement outcomes are possible as a result of the 
unpredictable introduction of different numbers of dark- 
counts. 
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BAYESIAN OPTIMISATION OF DETECTOR 
SIGNATURES 

We have argued (following Ref. [9j) that the operation 
of a number-resolving photo-detector can be fully char- 
acterized by the P{m\n) matrix. However, practically we 
are more interested in the converse conditional proba- 
bility, P{n\m). That is, what is the likelihood that the 
incident state has some number of photons, given a par- 
ticular measurement outcome. This in effect gives us the 
fidelity of the measurement result. 

This first step can be achieved using Bayes' Law to 
invert the conditional probability distribution. 



P{n = no\m = mo) = 



P{m = mo\n = no)P{n = no) 

P{m = mo) 
P{m. = mo|n = no)P{n = no) 
^^ P{m. = m.o\n = i)P{n = i) 

(3) 



Next we wish to optimise our signature. Specifically, 
for a given raw measurement signature m, we define the 
n which maximises P{n\m.) as the optimised signature. 
So, for example, if our detector gives us the measure- 
ment outcome r?i = 2, and then we find that P(n|2) is 



FIG. 1: (Colour online) P{m\n) matrices for a ideal, lossy 
and dark-county number-resolving detectors. ?7i = n = is 
front-left in each inset. 



maximised for n = 3, then we define mopt = 3 as the 
optimised signature. Thus, by reinterpreting m = 2 as 
^opt = 3, we have improved the confidence of our mea- 
surement outcome. Specifically, 



^opt(^) = argmax^P(n|m). 



(4) 



This idea was first employed in Ref. [7]. It can easily be 
seen that the state of highest confidence will always be a 
Fock state. 



EXAMPLE APPLICATIONS 

As an example test-bed for our technique, we consider 
two scenarios: (1) heralded preparation of Fock states 
using conditioned PDC sources, which is the staple for 
preparing heralded Fock states in present-day single pho- 
ton experiments, and (2) a uniform number distribution. 



The state emanating from a PDC can be expressed as 
a power series, 



IV'PDc) = v/r^^^x"|n).|n)i 



(5) 



where s and i are the signal and idler modes respectively. 
Thus, measuring the signal mode, the expected photon 
number distribution is given by 



P(n) = (1 - x')x 



2n 



(6) 



The photon number distribution for x = O-'^ is illus- 
trated in Fig. [2J For the uniform distribution we assume 
P{n) =0.1 in the range ... 9, otherwise 0. 
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FIG. 2: (Colour online) Photon number distribution, P(n), 
for a PDC state with x — 0-7- 

In Fig. Is] we plot the P{n\m) matrices, and in Fig. p] 
the corresponding optimisation maps (i.e. the mapping 
between the raw m and mopt) for our detector, assuming 
incident PDC and uniform states. When pioss = Pbc = 
the detector behaves ideally and the optimisation map 
is the trivial identity map, i.e. the raw detector signa- 
ture should be interpreted as is. As losses are introduced, 
Pioss > 0, it becomes desirable to compensate for the 
losses by incrementing the raw detector signature. On the 
other hand, when dark-counts are introduced, poc > 0, 
the opposite holds, and it becomes desirable to compen- 
sate for the dark-counts by decrementing the raw signa- 
ture. In the limit of large dark-count rates we observe 
that it is optimal to always assume that mopt = 0, since 
the detector is not giving us any useful information and 
the best we can do is guess that we have the most likely 
photon number, which in the case of our PDC source is 
n = 0. Note that the degree of required compensation in 
the measurement signature is strongly dependent on the 
underlying photon-number distribution. 

Importantly, in the regime where the optimisation map 
is non-trivial, the P{n\m) matrix is very 'washed-out'. 
For high fidelity measurements we desire highly localised 
peaks in the P{n\m) matrix. However, this does not occur 
in the regime where optimisation may take place. Thus, 
the applicability of the technique is rather limited. 

In Fig.|5]we compare the difference between the fidelity 
of the raw and optimised detector signatures with PDC 
and uniform input states. In the case of pioss = Pbc = 
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FIG. 3: (Colour online) P{n\m) matrices for ideal, lossy and 
dark-county detectors, and an incident PDC state. ??i = n = 
is front-left in each inset. 



we have ideal detector operation and there is no optimisa- 
tion. However, as losses and dark-counts are introduced, 
the optimisation technique improves the measurement fi- 
delity. The bias in the fidelity plots are very distinct for 
different photon-number distributions. 



CONCLUSION 

In conclusion, we have described a simple technique 
which may be employed to improve the confidence of 
single shot signatures generated from number-resolving 
photo-detectors. With sufficiently high loss or dark-count 
rates the technique can improve the confidence of the de- 
tection result. However, this implies that the technique 
only works in the regime where the detector is poor. 
This places limitations on when this approach might be 
employed. However, many present-day experiments oper- 
ate in the poor detector regime and this technique may 
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FIG. 4: Optimisation maps for ideal, lossy and dark-county 
detectors with PDC (left) and uniform (right) input states. 
Black boxes illustrate the mapping between raw and opti- 
mised signatures, while grey represents the identity (ideal) 
map and is used to guide the eye. m = mopt = is top-left in 
each inset. 
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FIG. 5: (Colour online) Raw (blue) and optimised (green) 
measurement fidelities for ideal, lossy and dark-county de- 
tectors with a PDC (left) and uniform (right) input number 
distributions. 



nonetheless find applicability. The example we presented 
was for the trivial case of an experiment consisting of 



a single photo-detector. However this technique can eas- 
ily be generalised to multiple detectors using the same 
Bayesian techniques. Additionally, any detector model 
can be employed whereby P{m\n) can be characterised, 
extending the applicability of this approach to many de- 
tector technologies, such as visible light photon counters 
(VLPCs) Uni E], time-multiplexed detectors (TMDs) 
J6l [T2Hl5] or other multiplexed detection approaches. We 
demonstrated that the optimisation map, which deter- 
mines how detector signatures should be compensated, 
is highly dependent on both the detector parameters and 
the underlying photon- number distribution. Thus, the 
applicability of this technique will be highly dependent 
on the specific application. 
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